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We investigate the solution of the continued fraction equation by which we determine 
"the renormalized angular momentum parameter", v, in the formalism developed by Leaver 
and Mano, Suzuki, and Takasugi. In this formalism, we describe the homogeneous solutions 
of the radial Teukolsky equation, which is the basic equation of the black hole perturbation 
formalism. We find that, contrary to the assumption made in previous works, the solution, 
. v, becomes complex valued as ui (the angular frequency) becomes large for each I and m 

(the degree and order of the spin- weighted spheroidal harmonics). We compare the power 
radiated by gravitational waves from a particle in a circular orbit in the equatorial plane 
around a Kerr black hole in two ways, one using the Mano- Suzuki- Takasugi formalism with 
complex v and the other using a direct numerical integration method. We find that the two 
»' | ' methods produce consistent results. These facts prove the validity of using complex solutions 

(3J[)! to determine the homogeneous solutions of the Teukolsky equation. 
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§1. Introduction 



Inspirals of stellar-mass compact objects into a supermassive black hole at galac- 
tic nuclei are expected to be one of the most important sources of gravitational 
waves for space-based detectors, such as the Laser Interferometer Space Antenna 
(LISA). 1 )' 2 ) Through the observation of gravitational waves from such systems, we 
may be able to obtain information regarding the central black hole's spacetime geom- 
etry encoded in multipole moments. We may also obtain astrophysical information, 
such as the mass distribution of compact objects in galactic nuclei. Theoretical 
waveforms of gravitational waves from such systems are required in order to extract 
such information through data analysis. 

To predict the waveforms of extreme mass ratio inspirals, we adopt the black hole 
perturbation approach. In this approach, a compact star is treated as a point particle, 
and the mass of the compact star, fj,, is assumed to be very small compared to the 
mass of the black hole, M, i.e. fj,/M <C 1. In this context, the Teukolsky equation 3 ) 
describes the evolution of a perturbation of the Kerr black hole spacetime. The 
standard approach to solve the Teukolsky equation is based on the Green function 
method. The Green function is constructed from two kinds of homogeneous solutions. 
The solution of the Teukolsky equation is obtained by the integration of the Green 
function multiplied by the source term. 

In a previous paper 4 ) (paper I), we discussed a numerical method to evaluate 
the homogeneous solutions using a formalism developed by Mano et al. 5 ) (MST). 
In this formalism, homogeneous solutions are expressed in series of hypergeometric 
functions around the horizon and in series of Coulomb wave functions around infinity. 
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This formalism is an extension of the formalism of Leaver, 6 ) in which only the series 
of Coulomb wave functions are given. One of the most important problems in the 
application of the Leaver and the MST formalisms to numerical computation is to 
determine the so-called "renormalized angular momentum", v, by solving an equa- 
tion of the form g{v) = 0, where g{v) is expressed in terms of continued fractions. 
We found that, for each s,£,m and q (where s is the spin index of the Teukolsky 
equation, t and m are the indices of the spin weighted spheroidal harmonics, and q 
is the Kerr parameter divided by the mass of the black hole, q = a/M), there is a 
maximum value of to for which we can obtain v, if we assume that the v is real. 

In this paper, we investigate in more detail the properties of the solution of the 
continued fraction equation g(v) = 0, defined in the Leaver and the MST formalisms. 
We search for v in the complex region if we cannot find v in the real region. We 
find that the solution is actually complex valued when no real solutions are found. 
Further, we find that the real part of the complex solutions, when they exist, always 
take integer or half-integer values. We confirm that these complex solutions of the 
continued fraction equation can be used to determine the homogeneous solutions of 
the Teukolsky equation. We compare the power radiated by gravitational waves from 
a particle in a circular orbit in the equatorial plane around a Kerr black hole com- 
puted by the MST formalism with complex v and by a direct numerical integration 
method. We find that the two methods give consistent results. These facts show the 
validity of using the complex solutions to determine the homogeneous solutions of 
the Teukolsky equation. By using complex z/, we can evaluate gravitational waves 
even in the case that Mui becomes large. This range of Muj is important for the 
LISA data analysis which deals with the case that a particle moves in the strong 
gravitational field or more generic orbits. 

This paper is organized as follows. In § 2, we briefly review the MST formalism 
and introduce the continued fraction equation, g(v) = 0. In § 3, we discuss the 
properties of the solutions of the continued fraction equation. In § 4, we compare 
the power radiated by gravitational waves from a particle around a Kerr black hole 
computed by the MST formalism and by a direct numerical integration method. 
Section 5 is devoted to a summary and discussion. Throughout this paper, we use 
units in which G = c = 1. 

§2. Analytic solutions of the homogeneous Teukolsky equation, and the 

continued fraction equation 

In this paper, we consider only the case oj > 0. Solutions in the case u < 
can be obtained from those in the case u > by the symmetry of the Teukolsky 
equation. 

The homogeneous Teukolsky equation is given by 
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where the potential term V(r) is given by 



K 2 + M(r — M)K , , , 

V(r) = y - — + 8iu;r + X. (2-2) 



Here A = r 2 - 2Mr + a 2 = (r - r + )(r - r_), with r± = M ± \/M 2 - a 2 and 
if = (r 2 + a 2 )u; — ma, and A is the eigenvalue of the angular Teukolsky equation. 

We consider the incoming solution of the homogeneous Teukolsky equation de- 
fined by the asymptotic behavior 

( Dtrans/\2 -ikr* f nr ~ . ~ , 

^•Jmw ^ r 3 B ^ e iur* + r -l j ginc^ e -^r* for f _^ +QO) I ^ 

where k = 00 — ma/2Mr + and r* is the tortoise coordinate defined by 

2Mr + , r — r_i_ 2Mr_ , r — r_ , 

r* = r-\ —In In . (24) 

r+ - r_ 2M r+ - r_ 2M V ; 

Here r± = M ± \J M 2 - a 2 . 

In the MST method, the homogeneous solutions of the Teukolsky equation are 
expressed in terms of two kinds of series of special functions. 5 )' 7 ) One consists of 
series of hyper geometric functions, and the other consists of series of Coulomb wave 
functions. The former is convergent at the horizon, while the latter at infinity. 

For the homogeneous solution, Rf^, we define p m by 

^in = e i, Kx ^_ x ys-i(e+r)/2^ _ x y(e~r ^ _ (3.5) 

The function p m is expressed in the form of a series of hypergeometric functions as 



Pi: 



= cinF(n + v + 1 — it, —n — v — it; 1 — s — ie — it; x), (2-6) 



where x = uj(r + — r)/en, e = 2Muj, k = \J\ — q 2 , q = a/M and r = (e — mq) / k, and 
F(a, (3; 7; x) is the hypergeometric function. 

Next, we present a solution in the form of a series of Coulomb wave functions. 
Let us denote a homogeneous solution by Rc- We define the function f v {z) through 
the relation 

^ = ^-.(l-^)-* +T|/2 ,„ W . (2-7) 
We express the function f u (z) in the form of a series of Coulomb wave functions as 

U(z) = £ i-ir tlllll^ ^n^i-is -e,z), (2-8) 

\V -\- 1 S ~r ie)n 

where z = tu(r — r_) and (a) n = r(a + n)/r(a), and Fn(t],z) is a Coulomb wave 
function defined by 

F N ( V , z) = e~ iz 2 N z N+1 F f^N + 2) ) ^ N + l ~^ 2N + 2 5 ^ ( 2 ' 9 ) 
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Here <P(a, (3; z) is the confluent hypergeometric function, which is regular at z = 
(see § 13 of Ref. 8)). The function R c is related to R in by 

R hl = K U R U C + K_ u _!R c v -\ (2-10) 

where K u is given by Eq. (2-20) of paper I. 

The expansion coefficients of the series of hypergeometric functions and the series 
of Coulomb wave functions, {a n }, satisfy the same three-term recurrence relation, 
given by 



<On+l + P>n + 7>n-l = 0, (2-11) 



where 



a v n = ieK(n + v)(2n + 2v - l)(n + v + 1 + s + ie) 

x(n + u + l + s- ie)(n + v + l + ir), (2-12) 
I3 y n = (-A - s(s + 1) + (n + i/)(n + 1/ + 1) + e 2 + e(e - mg)) 

x (n + i/)(n + 1/ + l)(2n + 2i/ + 3)(2n + 2i/ - 1) 
+e(e - m 9 )(s 2 + e 2 )(2n + 2^ + 3)(2n + 2u- 1), (2-13) 
7 £ = -ieK(n + u + l)(2n + 2^ + 3)(n + z/ - s + ie) 

x(n + — s — ie)(n + z/ — it). (2T4) 

The definitions of a„,/3^ and 7^ are different from those given in paper I by the 
factor (n + u)(n + v + l)(2n + 2zv + 3)(2n + 2i/ — 1). The series of Coulomb wave 
functions and the three term recurrence relation Eq. (2-11) are equivalent to those 
obtained by Leaver. 6 ) 

We note that the parameter v introduced in the above formulas does not exist in 
the Teukolsky equation. This parameter is introduced so that both series converge 
and actually represent a solution of the Teukolsky equation. Here, we review the 
method to determine v and obtain the convergence of both series expansions, based 
on Ref. 9). 

There are two independent solutions of Eq. (2-11) whose asymptotic behavior 
for n — > 00 are given by 

(1) • (2) o- 

hm n—— = — , hm — —- = — . (2-15) 

There are also two independent solutions whose asymptotic behavior for n — > — 00 
are given by 

(3) • (4) o- 
nm n , . = , lim — t-t— = . (2Tb) 

a n-l na n-l 



It is shown that in order to have convergence of Eqs. (2-6) and (2-8), the coefficients, 
dn \ must be used for n — > 00, and the coefficients must be used for n — > —00. 
These two solutions are called the "minimal solution" as n — >• ±00, respectively. 
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In fact, it can be proved that if we use these solutions as the expansion coefficients 
{o„}, the series of hypergeometric functions Eq. (2-6) converges for x in the range 
— oo < x < 0, and the series Eq. (2-8) converges for z > en, or equivalently, 
r > r+. However, these two solutions do not coincide in general. In order to obtain 
a consistent solution, the parameter u is used. 

We next introduce the following quantities: 

(1) (3) 

R n = —^-, L n = ^ . (2-17) 

a n-l a n+l 

We can express R n and L n in terms of continued fractions as 



Rn — 



in a n r n+ i a n+1 r n+2 



fin ftn+l @n4 



L r , 



J n+l f J n+2~ 



< <-l7n <-27n-l 
Pn~ Pn-1~ @n~2~ 



(2-18) 



(2-19) 



The convergence of the above continued fractions is guaranteed, because {a^} and 

{ai 3 '*} constitute minimal solutions. 10 ) 

In order to have the coincidence of the two solutions {a^} and {al 3 '}, we impose 
the relation 

RnL n -i = 1. (2-20) 

If we choose v such that this implicit equation is satisfied for a certain n, we can 
obtain a minimal solution {a n } that is valid over the entire range — oo < n < oo. 

In place of Eq. (2-20), we define an equivalent equation. First we divide Eq. 
(2-11) by a n and replace a n /a n ~i by R n and a n /a n +i by L n . We then find 

g n {u) = P» + c£R n+1 + YnLn-i = 0, (2-21) 

where R n +i and L n -\ ar e given by the continued fractions Eq. (2-18) and Eq. 
(2-19), respectively. Equation (2-21) is a transcendental equation expressed in terms 
of continued fractions. We call Eq. (2-21) the "continued fraction equation" in this 
paper. 



§3. Properties of the solutions of the continued fraction equation 

We first note that g n {y) i n Eq. (2-21) is invariant under the simultaneous trans- 
formation v <-> — v — 1 and n <-> —n. Thus, goiy) is symmetric under v <-> — v — 1. 
In particular, if we plot go{v) as a function of v, go(v) is symmetric about the line 
v = —1/2 in the complex plane. 
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Next, we note that the coefficients a^, fin and 7^ in Eqs. (2-12) - (2-14) contain 
v only in the form n + u. This implies that g n (v + l) = g n+ \{y). When v is a solution 
of Eq. (2-21), we have g n (v) = 0. This can be rewritten as 

9n+l{y) = fin+1 + a n+ iR n+2 + 7n+l^n = 0. 

We then have = g n+ \(y) = g n {y + 1). Thus, v + 1 is also a solution. Namely, if 
9n{ v ) = 0) 9n{y + k) = for an arbitrary integer k. These considerations show that 
it is sufficient to investigate g$(v) = to find the solution v. 

When e = 2Mw is small, there is an analytic expression of a solution v in the 
form of a series in e, given by 



v = I + ' 



21 + 1 
(I 2 - s 2 ) 2 
' (21 - 1)21(21 + 



2 * 2 , [(^ + 1) 2 -* 



1(1 + 1) + (2Z + l)(2Z + 2)(2Z + 3) 

e 2 + 0(e 3 ). (3-1) 



Thus, when e — > 0, takes an integer value u = I + k (where k is an arbitrary 
integer). As an example, we consider the case I = m = 2, q = and k = 0. Then, 
as w increases from 0, the solution z/ decreases from I and approaches 1 — 1/2 near 
Mw = 0.36. If Mw is larger than 0.36, the real solution v disappears. The plot 
of go(x)/l 5 is shown in Fig. 1 in the case Muj < 0.36 and in Fig. 2 in the case 
Mlj > 0.36. The factor Z™ 5 is introduced so that the function go(x)/l 5 remains order 
unity. At first glance, in Fig. 2, there seems to be a solution at x = I — 1/2 = 3/2, 
x = 1 or x = 2. However, at these points, go(x) always becomes 0, irrespective of 
ll>. These are solutions of go(x) = mathematically, but they are not solutions that 
can be used to connect two minimal solutions in the limits n = ±oo, discussed in 
the previous section. In Table I of paper I, we listed the maximum values of Muj for 
which real v is found. We also list them in Table I for convenience. We express the 
maximum value of Muj as (Mw) max . 

Next we consider the solution of go(z) = in the complex z plane in the case 
that Muj is larger than (Mw) max and we cannot find a real solution. In Fig. 3, we 
plot a contour of |<7o(z)| i n the case Muj = 0.5. We find that there is a minimum near 
Re(z) = 3/2, lm(z) = 0.36. In Fig. 4, we plot Re(go(z)) and lm(go(z)) as functions 
of lm(z) at Re(z) = 3/2. It is evident that there is a solution near z = 3/2 + 0.36i. 
The precise value of the solution is z = 3/2 + 0.36188061539416i. It is also suggested 
by Fig. 3 that there are no solutions other than this value. We find that when there 
are no real solutions, a complex solution always exists. Further, the real part of the 
complex solution is always a half-integer or integer. These properties are the same 
when the parameters, s,l,m and q have different values. 

We consider the case in which Muj is much larger than (i\/fw) max . In Table II, 
we list numerical values of the solution v in the case s = —2, I = m = 2 and q = 0. 
We find that Im(z/) increases as Muj increases, and it reaches a maximum value near 
Muj = 0.53. Then, Im(zv) begins to decrease and approaches 0, and the real solution 
appears again at Muj = 0.594. This new real solution starts at the half integer value 
v = 3/2, and it approaches the integer value v = 2. At Muj = 0.604, where the 
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Table I. Maximum values of Muj for which real solutions v are found. 
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real solution reaches v = 2, the real solution disappears, and the imaginary part 
appears again. In Table III, we list numerical values of the solution v in the case 
s = —2, I = m = 2 and q = ±0.99. When q = 0.99, the real solution disappears 
at Muj = 0.39, and a complex solution with Re(v) = 3/2 appears. Then lm(u) 
continues to increase until Mu> = 3.0. When q = —0.99, the real solution disappears 
at Muj = 0.33, and a complex solution with Re(^) = 3/2 appears. Then Im(i/) 
continues to increase until Muj = 3.0. In Tables II and III, we list values for the case 
Muj < 3.0. Above this value, it becomes difficult to evaluate the continued fraction 
equation numerically, because of the loss of precision. This happens because as Mui 
becomes large, (3 n approaches — (a^i? n +i +7^L n _i) irrespective of v. 

In Figs. 5-13, we plot the real and imaginary parts of v for I = m = 2, 3, 4 
and q = 0, ±0.99. In Figs. 14 and 15, we plot the real and imaginary parts of v 
for I = 2, . . . , 9 in the case q = 0. We find that the properties of v are very similar 
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1 1.2 1.4 1.6 1.8 2 

V 



Fig. 1. The plot of g {x)/l 5 for Muj = 0.20 
and 0.36. Here s = —2, I = m = 2 and 
q = 0. There exist real solutions in these 
cases. 



4,5 




1 1.2 1.4 1.6 1.8 2 

V 

Fig. 2. The plot of go{x)/l 5 for Muj = 0.37, 
0.40 and 0.50. Here s = -2, I = m = 2 
and q = 0. No real solutions are found in 
these cases. 



, (v) 1/1 



1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 

Re(v) 

Fig. 3. The contour plot of \go(z)\/l in the 
case Muj = 0.50, s — —2, I — m = 2 and 
g = 0. There is a complex solution at v — 
1.5 +0.36188061539416i. 




/ \ : Relg (!))/! 

\ 

* 'Imlg (z))/l 

X 



Im(z) 

Fig. 4. The real and imaginary parts of the 
function go(z)/l at Re(z) = 1.5 in the 
case Muj = 0.50, s = —2, I — m — 
2 and q = 0. The real and imaginary 
parts become simultaneously at Im(z) = 
0.36188061539416. Thus, this function sat- 
isfies go(z) — 0. 



to those in the case 1 = 2, although there are several quantitative differences. For 
example, in the case I = 5, in Figs. 14 and 15, Re(^) reaches 4.5 at Mu = 0.776, 
and a complex solution appears. However, this complex solution exists only in a 
very small interval of Mu, 0.776 < Muj < 0.778, and a real solution appears again 
for Muj > 0.778. 

In Tables IV- VI, we list the values of v for various Muj in the cases s = 0, 
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I = m = 2 and q = 0, ±0.99. We find behavior of v very similar to that in the case 
s = -2. 

Here, we briefly explain our numerical method of searching for a solution of 
go(z) = 0. When we search for a real solution, it should be sufficient, in principle, to 
search the region I — 1/2 < x < I, because of the symmetry property of the solution 
mentioned above. However, in some cases, the function go(x) becomes very steep 
around the solution in that region. In such cases, it is often easier to search for the 
solution in the region I — 1 < x < I — 1/2. We thus search the region I — 1 < x < I. 
We discard solutions at integer and half-integer values, x = I — 1,1 — 1/2,1. 

When we do not find any solutions on the real z axis, we search for a solution in 
the imaginary direction. We consider only half-integer and integer values of Re(z), 
and search in the imaginary direction of the z-plane. Thus, the problem becomes a 
simple one-dimensional search for a root of the complex equation go(z) = 0. Further, 
we utilize a very useful property of go(z): The function go(z) (with complex z) is 
real when z = — 1/2 + iy (with real variable y). This is derived using the fact that 
Pk = P^k an d «fc7£i+i = ipi-k-il-k)* (where * denotes complex conjugation, and 
k is an arbitrary integer). Thus, when we search for a solution with half-integer 
real part, it is convenient to search at Re(z) = —1/2, since when v is a solution, 
v + k (with k an arbitrary integer) is also a solution. For a solution with integer 
real part, we search for a solution which simultaneously satisfies Ke(go(z)) = and 
lm(go(z)) = 0. In particular, we search for the real part Re(z) = 1 irrespective of 
I, because the function go{z) becomes steeper around the solution as Re(z) becomes 
larger, and it becomes more difficult to find the solution in such cases. 

§4. Gravitational wave luminosity 

Since the role of v is to match two minimal solutions in the limits n — ► ±oo of 
the recurrence relation, it is evident that once we find a solution of Eq. (2-21), it 
can be used to generate homogeneous solutions irrespective of whether it is real or 
complex. In this section, we confirm this by using the luminosity of gravitational 
waves radiated by a particle in a circular, equatorial orbit around a Kerr black hole. 
The complete formulas are given in Appendix A of paper I. We compare the numerical 
data with those computed with the numerical integration method of Kennefick. *) 

The innermost stable circular orbit (ISCO) around a Kerr black hole is given by 

nsco = M[3 + Z 2 T y / {3-Z 1 ){3 + Z l + 2Z 2 )}, (4-1) 
Zi = 1 + (1 - g 2 ) 1/3 [(l + g) 1/3 + (1 - q) 1/3 }, (4-2) 

Z 2 = yj3q* + Zl (4-3) 
where the upper and lower signs refer to the cases q > and q < 0, respectively. 



*' The numerical data were kindly provided to us by D. Kennefick and are based on the work 
presented in Ref. 11). 
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Re(v) 
Im(v) 



Re(v) 
Im(v) 



Fig. 5. The values of the real and imaginary 
parts of v as functions of Muj in the case 
s — —2, I = m = 2 and q — 0. 



Fig. 6. The values of the real and imaginary 
parts of v as functions of Muj in the case 
s = -2, I = m = 2 and q = 0.99. 




Fig. 7. The values of the real and imaginary parts of v as functions of Muj in the case s = —2, 
I = m = 2 and q = -0.99. 



The orbital angular frequency of a circular orbit is given by 



M 1 ' 2 



r 3 / 2 (l + g(M/r) 3 / 2 )' 



(4-4) 



The angular frequency of gravitational waves corresponding to the harmonic m is 
given by oj = mfl. For q > 0, the orbital radius of the ISCO becomes smaller than 
that in the Schwarzschild case, and it approaches the event horizon as q — > 1. Thus, 
Muj around the ISCO can be larger than (Mw) max . 

In Table VII, we list the numerical data for the gravitational wave luminosity 
induced by a particle in a circular orbit in the case that Muj is larger than (Mw) max . 
We find that the numerical data computed using a complex v agree with those 
computed using the numerical integration method to about 5 or 6 digits. This 
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Fig. 8. The values of the real and imaginary 
parts of v as functions of Muj in the case 
s — —2, I = m = 3 and q — 0. 



Fig. 9. The values of the real and imaginary 
parts of v as functions of Muj in the case 
s = -2, I = m = 3 and q = 0.99. 




Fig. 10. The values of the real and imaginary parts of v as functions of Mlu in the case s 
I — m — 3 and q = —0.99. 



is approximately the same degree of accuracy of the numerical integration. This 
demonstrates the validity of using a complex v to produce the homogeneous solutions 
of the Teukolsky equation. 

In the case of eccentric orbits, there are infinite number of higher harmonics for 
each l,m. Thus, in principle, there are many harmonics for which Muj is larger than 
(Mw) max . Here, we estimate the range of harmonics needed to express the total 
power of the gravitational waves with the relative error less than 10 -5 using the 
Newtonian, quadrupole approximation, in the case of a non-spinning binary. The 
power radiated in nth harmonics is given by 12 ) 

32 m\rn\[rn x + m 2 ) 
P{n) = — g g(n,e), (4-5) 

r a 
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Fig. 13. The values of the real and imaginary parts of v as functions of Mlo in the case s = —2, 
I = m = 4 and q = -0.99. 



where mi and mi are the masses of the two stars forming the binary, e is the 
eccentricity, r a is the semi-major axis of the orbit, and we have 

2 I 2 

J n _ 2 (ne) - 2eJ n _i(ne) + -J n (ne) + 2eJ n+ i(ne) - J n+2 (ne) 

re 

+ (1 - e 2 )[J n - 2 (ne) - 2J n (ne) + J„ +2 (ne)] 2 + ^( J nM) 2 j> C 4 " 6 ) 

where J n (z) is the Bessel function. The angular frequency for each harmonic is given 
by oj = nujo, with u>o = (mi + mz) 1 / 2 /^ ■ In Fig. 16, we plot the contours of the 
values of Mlo that must be included in the total power as functions of e and r a . We 
find that for small r a , there is a relatively wide range of values of e for which the 



g(n,e) 
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Real part o§ 




M to 

Fig. 14. The values of the real parts of v as 
functions of Mu in the cases I = 2, 3, • • ■ ,9, 
s = —2 and q = 0. 



Imaginary part ctf 







1=3 




1=4 




1=5 

. 1=7 


s * 

/ J'/ 


1=8 




1=9 


/ ' ' '* / \ .' 
















s / '' /t ■■■■ ' - < ' 




/ ' ■"■«* / r ' 1 




/ .'</ V / ;' 

/ • • ,,// ....I , 




' ""J ' ni 'a > 



M to 



Fig. 15. The values of the imaginary parts of 
v as functions of Mui in the cases from I = 
2,3, ••• ,9, s = -2 and q = 0. 




B D.l D.2 D.3 D. 4 D.5 0.6 0.7 0.0 0.5 

eccentricity 



Fig. 16. The contours of the values of Mu){= Mnu>o) that must be included in the evaluation of 
the total power in order to obtain the the relative error less than 10 -5 . Here, M = mi + m2. 
The thick curve corresponds to the contour for Muo — 0.36, which is equal to (Mtj) max , in the 
case I = m = 2 and q — 0, above which v becomes complex. The region under thick dashed 
curve should be excluded from consideration, because there, the pericenter distance is less than 
2M. Thus, in the region between the thick curve and the thick dashed curve, the parameter v 
can have a complex value. 



required value of Moj exceeds (Mw) max = 0.36 in the case I = m = 2 and q = 0. 
In such cases, complex v must be employed in order to evaluate the gravitational 
radiation. Of course, we should keep in mind that, because this is based on the 
Newtonian, quadrupole approximation, this result will be changed quantitatively in 
the case of fully relativistic treatment. 
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§5. Summary and discussion 

We have investigated the numerical properties of the solutions of the continued 
fraction equation derived by Leaver and Mano et al. in the formalism of the analytic 
representation of the homogeneous solution of the Teukolsky equation. We have 
found that, for each s,l,m and q, the solution v takes a real or complex value, 
depending on the value of to. We have also found that, when v is complex, its real 
part is always an integer or half integer. Although there are certain regularities in the 
behavior of the solution, we have not determined how such behavior is controlled by 
the continued fraction equation. We hope to investigate the behavior of v analytically 
in the future. 

In order to confirm the existence of complex u, we have computed the gravita- 
tional wave luminosity emitted by a particle in circular, equatorial orbit around a 
Kerr black hole in the case that v is complex, and compared the results with those 
computed using a numerical integration method. We found that the two results are 
consistent. This shows clearly the validity of the existence of complex v. Even if 
Muj becomes larger than a certain critical value and a real v does not exist, we can 
use the MST formalism by treating v as a complex value. Thus, the region of oj in 
which the MST formalism can be used becomes much wider than in the case that 
we consider only real u, as in paper I. 

There is, however, a limitation on the values of lj for which we can use the 
MST formalism for actual numerical computations. When Muj becomes very large 
(say Muj > 5 in the case s = —2, I = 2 and q = 0), it becomes very difficult 
to evaluate the continued fraction equation (2-21) accurately because of the loss of 
accuracy. As oj becomes large, (3 n approaches — {a^Rn+i + Tn-^n-i) irrespective of 
v. Thus, it becomes difficult to evaluate Eq. (2-21) numerically with high accuracy. 
We must solve this problem because Moj becomes very large in the cases of generic 
orbits, which are important for the data analysis for LISA. A possible solution is to 
derive the high frequency expansion of v. If we cannot solve the problem, we may 
need to use numerical integration methods to evaluate the homogeneous solutions 
of the Teukolsky equation when Muj becomes very large. However, we believe that 
the MST method supplemented by some numerical integration is a very powerful 
method for the numerical computation of the Teukolsky equation. 

Based on the analysis presented in this paper, we plan to apply the MST formal- 
ism to the evaluation of the gravitational wave luminosity in the cases of eccentric 
orbits and more generic, non-equatorial orbits, and to investigate the data analysis 
issue for LISA in the near future. 
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Table II. Table of v for various values of Mui in the case s = —2, I — m — 2 and q = 0. 



Muj 


Re(z/) 


Im(i^) 





1000 


1 


9793154547208 





0000000000000 





2000 


1 


9129832302687 





0000000000000 





3000 


1 


7792805424199 





0000000000000 





3100 


1 


7594378913945 





0000000000000 





3200 


1 


7375329132632 





0000000000000 





3300 


1 


7129796083608 





0000000000000 





3400 


1 


6847877144372 





0000000000000 





3500 


1 


6510165965479 





0000000000000 





3600 


1 


6066167352463 





0000000000000 





3700 


1 


5000000000000 





0117805002266 





3800 


1 


5000000000000 





1080591170505 





4000 


1 


5000000000000 





1862468531447 





4200 


1 


5000000000000 





2393542510924 





4400 


1 


5000000000000 





2810676246093 





4600 


1 


5000000000000 





3149109523037 





4800 


1 


5000000000000 





3419264795996 





5000 


1 


5000000000000 





3618806153941 





5200 


1 


5000000000000 





3732918843902 





5300 


1 


5000000000000 





3748108599871 





5400 


1 


5000000000000 





3725753632958 





5600 


1 


5000000000000 





3508153681782 





5800 


1 


5000000000000 





2767032401575 





5920 


1 


5000000000000 





1112238120363 





5930 


1 


5000000000000 





0654237590662 





5940 


1 


5671748369160 





0000000000000 





5960 


1 


6572267840117 





0000000000000 





5980 


1 


7231093132436 





0000000000000 





6000 


1 


7878302655744 





0000000000000 





6020 


1 


8649828581648 





0000000000000 





6030 


1 


9221344236998 





0000000000000 





6040 


2 


0000000000000 





0739848339575 





6060 


2 


0000000000000 





1644030708573 





6080 


2 


0000000000000 





2167770350574 





7000 


2 


0000000000000 





8003377636925 





8000 


2 


0000000000000 


1 


1099466644118 





9000 


2 


0000000000000 


1 


3699138540831 


1 


0000 


2 


0000000000000 


1 


6085538776570 


1 


2000 


2 


0000000000000 


2 


0541150948294 


1 


4000 


2 


0000000000000 


2 


4772658891494 


1 


6000 


2 


0000000000000 


2 


8879097347262 


1 


8000 


2 


0000000000000 


3 


2903727081665 


2 


0000 


2 


0000000000000 


3 


6867890278893 


2 


2000 


2 


0000000000000 


4 


0782125805902 


2 


4000 


2 


0000000000000 


4 


4650426712959 


2 


6000 


2 


0000000000000 


4 


8471678988185 


2 


8000 


2 


0000000000000 


5 


2239336390821 


3 


0000 


2 


0000000000000 


5 


5939000509184 



New Numerical Methods to Evaluate the Homogeneous Teukolsky Equation 



Table III. Table of v for various values of Mlu in the case q — 0.99 and —0.99, s — —2, I = m 







q = 


0.99 




q = 


-0.99 


Muj 


Re(i/) 


Im(i/) 


Re(z/) 


Im(z/j 


n 
u 


i o 


i 
i 


o. 80£ 9 ^ ^ ^008 qq 


u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


Q778£81 QOQ/1 1 Q 

y/ / ouoiyuy^iy 


u 


0000000000000 

uuuuuuuuuuuuu 


u 


zu 


i 
i 


09^0^1 8*3/1/10*39 


u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


808/19^77^978*3 

oyo^zo i i oz / oo 


n 
u 


0000000000000 

uuuuuuuuuuuuu 


u 


*30 


i 
i 


81 8*3^81 8*3*3078 
olooOolooou f o 


u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


£8^791 88£007Q 
DoO / ZlooDUU i y 


n 
u 


0000000000000 

uuuuuuuuuuuuu 


n 
u 


*31 


i 


8097 c if : i8870f : i1 ^ 
OUZ ( uUOO / uu±o 


o 

u 


0000000000000 
uuuuuuuuuuuuu 


i 
i 


uoi loozouyuoz 


n 

u 


0000000000000 
uuuuuuuuuuuuu 


n 
u 


*39 

oz 


i 
i 


78^/1 £09088/1 ^/l 

i oo z J:Dyzyoo z ±0'i 


n 
u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


^/IQQQI £8770^0 

o^yyyiDo/ *yoy 


n 
u 


0000000000000 

uuuuuuuuuuuuu 


u 


*3*3 

JO 


i 
i 


7££0 c i1 9Q/I71 OK 

(uouoizy^i/ iuo 


n 
u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


1 901 A £0888^*31 
1ZU14DU000001 


n 
u 


04: 


i 


7/1*38/18/19 00 £0*3 


o 

u 


0000000000000 
uuuuuuuuuuuuu 


i 
i 


^nnnooooooooo 

ouuuuuuuuuuuu 


o 

u 


1 7880 c i7'37708/l 
± f ooyo foil U04: 


n 
u 


*3^ 


i 
i 


71 789/1 77Q1 9*3*3 
/ 1 loZ4( /yizoo 


n 
u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


99/108/100^1 8£9 
zz^ucfiuyo looz 


n 
U 


*3£ 


i 
i 


£8£1 7190^^9800 

ooonzuoozoyy 


u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


9£981 ^80^*379^ 
ZDZolOoUDo ( ZO 


n 
u 


*37 

O ( 


i 
i 


£/1/17Q/191 /I01 70 


n 
u 


0000000000000 

uuuuuuuuuuuuu 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


907^*3/1 *31 70*307 

zy * Do4oi / yoy / 


n 
u 


*38 
oo 


i 

-L 


^ 7 Q *3 Q £ £ ^ ^ ^ *3 
o i youuyuuoooo 


o 

u 


0000000000000 
uuuuuuuuuuuuu 


i 
i 


^nnnooooooooo 

ouuuuuuuuuuuu 


n 

u 


*39 0709 9/170*3/17 
ozy t yzz'i i uoi i 


n 
u 


*3Q 


1 
1 


^OOOOOOOOOOOO 

ouuuuuuuuuuuu 


u 


0Q7/I A 8Q0Q/I A *30 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


*3£001 *391 090/1^ 

ODUUloZlUZU z ±U 


u 


/in 

4U 


1 
1 


^000000000000 

ouuuuuuuuuuuu 


A 

u 


1 £981 ^90^9/17*3 
-Luzoiozyoz^i / o 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


'3887^8/11 *3/18^7 
OOO / U0410400 i 


n 
u 


^0 
OU 


1 
1 


^000000000000 

ouuuuuuuuuuuu 


u 


^*3£8/1Q81 £90/1 1 

ooOo^yoiozu^i 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


£*3£ c ;'3( : iOQ/IQQ99 

DoDOoDuy^yy zz 


n 
U 


£0 

OU 


1 
1 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


8£790£79£7Q7^ 
00 J ZUO / ZD i y J 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


n 
u 


8^71 9981 1 990/1 
OO ( lZZollZZU z ± 


n 
u 


70 
1 u 


1 
1 


p^oooooooooooo 
ouuuuuuuuuuuu 


1 


1 8*38£Pi1 OQO*370 
loooDOlUyUo / u 


i 


p;oooooooooooo 
ouuuuuuuuuuuu 


1 
1 


070/1 ££1 89^Q£/I 
U 1 U4DDloZ0yD4 


n 
u 


80 

oU 


1 
1 


^000000000000 

ouuuuuuuuuuuu 


1 

1 


/1 8/1 09££/1 8Q87£ 

lo^uzuoioyo / o 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


1 

1 


9891 9/1^1 ^^fiQQ 
zoz izio loooyy 


n 
u 


00 

yu 


1 
1 


^000000000000 

ouuuuuuuuuuuu 


1 

1 


7701 *38£ p i£'3708 
/ /UlooOOOoiUo 


i 

j. 


^000000000000 

ouuuuuuuuuuuu 


1 

1 


A Q *3 8 8 A Q *3 Q 79 A 7 

4yooo4yoy / Z4 ( 


± 


uu 


1 
1 


p^oooooooooooo 
ouuuuuuuuuuuu 


z 


O/l KA *3*30£9Q£99 
U404ooUDZyDZZ 


i 


p;oooooooooooo 
ouuuuuuuuuuuu 


1 
1 


70£1 *371 ^^9.9.(^1 
(UDloi lOoooDi 


i 
i 


1 o 

1U 


1 
1 


^000000000000 

ouuuuuuuuuuuu 


o 
z 


*31 9^/101 79Q£1 Q 
oizo^ui / zyoiy 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


1 

1 


01 879789881 £9 
y lo / Z / oZoolDZ 


i 
i 


90 

zu 


1 
1 


^000000000000 

ouuuuuuuuuuuu 


o 
z 


^7*3/1090*38008^ 


i 
i 


^000000000000 

ouuuuuuuuuuuu 


o 
z 


1*31*3*31 7fi77£98 
lolOOl (D( (DZo 


1 


30 


1 


5000000000000 


2 


8294329783970 


1 


5000000000000 


2 


3436142271920 


1 


40 


1 


5000000000000 


3 


0816629547228 


1 


5000000000000 


2 


5552922279595 


1 


50 


1 


5000000000000 


3 


3308588672551 


1 


5000000000000 


2 


7661482256677 


1 


60 


1 


5000000000000 


3 


5775977157128 


1 


5000000000000 


2 


9760238443012 


1 


70 


1 


5000000000000 


3 


8223206679586 


1 


5000000000000 


3 


1848072932141 


1 


80 


1 


5000000000000 


4 


0653697494531 


1 


5000000000000 


3 


3924204618486 


1 


90 


1 


5000000000000 


4 


3070135210021 


1 


5000000000000 


3 


5988078524059 


2 


00 


1 


5000000000000 


4 


5474653357747 


1 


5000000000000 


3 


8039277638612 


2 


10 


1 


5000000000000 


4 


7868965237373 


1 


5000000000000 


4 


0077454471779 


2 


20 


1 


5000000000000 


5 


0254460525305 


1 


5000000000000 


4 


2102277404164 


2 


30 


1 


5000000000000 


5 


2632277012574 


1 


5000000000000 


4 


4113386736802 


2 


40 


1 


5000000000000 


5 


5003354519402 


1 


5000000000000 


4 


6110355696752 


2 


50 


1 


5000000000000 


5 


7368475857391 


1 


5000000000000 


4 


8092651925867 


2 


60 


1 


5000000000000 


5 


9728298259083 


1 


5000000000000 


5 


0059594834575 


2 


70 


1 


5000000000000 


6 


2083377711225 


1 


5000000000000 


5 


2010303381768 


2 


80 


1 


5000000000000 


6 


4434187950814 


1 


5000000000000 


5 


3943626971054 


2 


90 


1 


5000000000000 


6 


6781135409542 


1 


5000000000000 


5 


5858048534288 


3 


00 


1 


5000000000000 


6 


9124571056503 


1 


5000000000000 


5 


7751542098469 
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Table IV. Table of v for various values of Mu; in the case s = 0, I — m — 2 and q = 0. 



Muj 


Re(u) 


Im(i^) 


n 
u 


i nno 

1UUU 


i 

j. 


yo^oz / oui^ooz 


n 
u 


nnnnnnnnnnnnn 
uuuuuuuuuuuuu 


n 
u 


9000 
zuuu 


i 
i 


Q^{7fi7489zL1 ^81 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


^ooo 

ouuu 


i 

j. 


8 c .999nfi W^fi ^ ^ 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


4000 

4UUU 


i 
i 


7nfiQ9' : !zl9ni 119 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


u 


/i i oo 


i 
i 


RSfil 87fi8fiQ' : !8n 
uouio i uouyoou 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 




i 
i 


fifi^l QzL^zLSfini 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


4^00 


i 
i 


fi'^fiQfiflS'^SI 8Q 
DoOyOOojioloy 


n 
u 


nnnnnnnnnnnnn 
uuuuuuuuuuuuu 


u 


/i /inn 


i 
i 


UUOo^Zuluy^ol 


n 
u 


nnnnnnnnnnnnn 
uuuuuuuuuuuuu 


u 


4^00 


i 

j. 


c .fins99n^ c i8fii r 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


4^00 


i 

j. 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


n ^ S n A 8 P\ R 9 A 4 4 ^ 


u 


/1700 


i 

j. 


nnnnnnnnnnnnn 
ouuuuuuuuuuuu 


u 


iuuiou z ±zoi z ±yu 


n 
u 


4800 


i 

j. 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


1 97^^' : i8f : 17Sl 9n 
±Z ( oooou / o±zu 


n 
u 


4qoo 


i 
i 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


i-*± i uozoy < zouo 


n 

u 


^ooo 

ouuu 


i 
i 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


1 ^^^9879^^894 


n 

u 


no 

OIUU 


i 

j. 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


1 7^1 8Qn77fin8f : i 
i i oioyu/ /uuou 


n 
u 


^900 
ozuu 


i 
i 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


1 8^47^947^7 


n 
u 


^oo 


i 
i 


nnnnnnnnnnnnn 
ouuuuuuuuuuuu 


u 


1 88n^1 R71 4fi71 


n 
u 


^400 

04:UU 


i 
i 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


1 88^Q8f : 17 c i4n7n 
loooyou / o^u / u 


n 
u 


^oo 

oouu 


i 
i 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


1 8^4^8^9981 ^9 


n 
u 


ouuu 


i 
i 


nnnnnnnnnnnnn 
ouuuuuuuuuuuu 


n 
u 


1 71 99^1 ^^7^8/1 
1 ( IZZOlOO ( OOT 


n 
u 


^700 
O I uu 


i 
i 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


1 4791 1 ^nn^i 87 

ill Zl lOUUOlO / 


n 
u 


^800 

oouu 


i 
i 


nnnnnnnnnnnnn 

ouuuuuuuuuuuu 


n 
u 


nQ74n^l 77^1 7^ 
uy /luoi / ( o± < o 


u 


^ooo 
oyuu 


i 

j. 


^Qfi79£fin* ; t7/in p ; 
oyu / zdouo / tuo 


U 


nnnnnnnnnnnnn 
uuuuuuuuuuuuu 


n 
u 


fiooo 

uuuu 


i 
i 


uy / oooyy±i i uu 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


8040 


i 
i 


7^70771 ^n^7nQ 
i o i u t 1 rooo / uy 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


U 


fi080 

UUoU 


i 
i 


77Qfi/l8ni *}£99n 

/ / yu^ooioDzzo 


n 
u 


nnnnnnnnnnnnn 
uuuuuuuuuuuuu 


n 
u 


fil 90 


i 
i 


R9Qni Q901 f)899 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


fii ^o 

UlUU 


i 
i 


RQ^W^nQSQQnnn 
oyoiOjyoyyouu 


n 
u 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


U 


^900 
ozuu 


9 
Z 


uuuuuuuuuuuuu 


n 
u 


n881 9^9889n97 
UOOlZOZOOZUZ ( 


u 


uuuu 


9 
Z 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


41 897n8m 87^^ 

'rtiOZ ( UOUIO / uo 


n 
u 


7000 
I uuu 


9 
Z 


nnnnnnnnnnnnn 

uuuuuuuuuuuuu 


n 
u 


c i7Q8 c iQn4Ql n9Q 


u 


8000 
oUUU 


9 
z 


nnnnnnnnnnnnn 
uuuuuuuuuuuuu 


n 
u 


oDozDtjzDlUyUU 


o 


9000 


2 


0000000000000 


1 


0958382861914 


1 


0000 


2 


0000000000000 


1 


2916407295318 


1 


1000 


2 


0000000000000 


1 


4596540676691 


1 


2000 


2 


0000000000000 


1 


5871735288388 


1 


3000 


2 


0000000000000 


1 


5757480036493 


1 


4000 


1 


5000000000000 


1 


8980859589471 


1 


6000 


1 


5000000000000 


2 


5246610689374 


1 


8000 


1 


5000000000000 


3 


0134306373653 


2 


0000 


1 


5000000000000 


3 


4652896348154 


2 


2000 


1 


5000000000000 


3 


8988668745912 


2 


4000 


1 


5000000000000 


4 


3212911554702 


2 


6000 


1 


5000000000000 


4 


7360933390574 


2 


8000 


1 


5000000000000 


5 


1452723630699 


3 


0000 


1 


5000000000000 


5 


5500484659006 



New Numerical Methods to Evaluate the Homogeneous Teukolsky Equation 



Table V. Table of v for various values of Mlo in the case s — 0, / = m = 2 and q = 0.99. 



Mlo 


Re(z/) 


Im(V) 





10 


1 


9855006468982 





0000000000000 





20 


1 


9433423104280 





0000000000000 





30 


1 


8738919793067 





0000000000000 





40 


1 


7745498379081 





0000000000000 





50 


1 


6249603205669 





0000000000000 





51 


1 


6016482354182 





0000000000000 





52 


1 


5720774742602 





0000000000000 





53 


1 


5140777829215 





0000000000000 





54 


1 


5000000000000 





0683318312500 





55 


1 


5000000000000 





0970789832967 





56 


1 


5000000000000 





1187106486632 





57 


1 


5000000000000 





1367752823053 





58 


1 


5000000000000 





1526950899429 





59 


1 


5000000000000 





1672593963838 





60 


1 


5000000000000 





1809982771760 





70 


1 


5000000000000 





3387878153683 





80 


1 


5000000000000 





5844007364833 





90 


1 


5000000000000 





8212346456388 


1 


00 


1 


5000000000000 





9679578324407 


1 


10 


2 


0000000000000 


1 


0719514201554 


1 


20 


2 


0000000000000 


1 


5606928653937 


1 


30 


2 


0000000000000 


1 


8822793013188 


1 


40 


2 


0000000000000 


2 


1581224329993 


1 


50 


2 


0000000000000 


2 


4081352814494 


1 


60 


2 


0000000000000 


2 


6392239009791 


1 


70 


2 


0000000000000 
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8536938783809 


1 


80 


2 


0000000000000 


3 
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1 


90 


2 


0000000000000 


3 


2246989210551 


2 


00 


2 


0000000000000 


3 


3561709666536 


2 


10 


2 


0000000000000 


3 


2904400721719 
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20 


1 
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3 
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30 


1 
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0811236068986 
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1 


5000000000000 
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2 


60 


1 


5000000000000 
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1 


5000000000000 
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0989848264810 
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1 


5000000000000 


5 


3251065286991 
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1 
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5 


5449488222337 


3 


00 


1 


5000000000000 


5 


7592716343452 
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Table VI. Table of v for various values of Mcu in the case s — 0, I = m = 2 and q — —0.99. 



Mlj 


Re(z/) 


Im(i') 





10 


1 


9840980928918 





0000000000000 





20 


1 


9308382464262 





0000000000000 





30 


1 


8203733727014 





0000000000000 





31 


1 


8041051082066 





0000000000000 





32 


1 


7863405245509 





0000000000000 





33 


1 


7667869072726 





0000000000000 





34 


1 
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0000000000000 
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1 
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1 
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1 
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1 
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1 


5000000000000 
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00 
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5000000000000 
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5000000000000 
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5000000000000 
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8121617032199 
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70 


1 


5000000000000 


5 


0301522240344 


2 


80 


1 


5000000000000 


5 


2450395718077 


2 


90 


1 


5000000000000 


5 


4571399180362 


3 


00 


1 


5000000000000 


5 


6666774973656 
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Table VII. Comparison of the energy fluxes obtained using the numerical integration method of 
Kennefick and the method presented here in the case ro = 1.55M and q = 0.99, with complex 
v. 



1 


m 




Im(!/) 


Numerical integ 


ration 


This paper 


Relative error 


2 


2 


1.5 


1 


1374192131794 


3 


568050135 


X 


10" 


2 


3 


568033154338851 


X 


10" 


2 


4.76 


X 


10" 


6 


3 


3 


2.5 


1 


4249576682707 


2 


152962111 


X 


10" 


2 


2 


152959342790158 


X 


10" 


-2 


1.29 


X 


10" 


6 


4 


4 


3.5 


1 


7677955367662 


1 


230541176 


X 


10" 


2 


1 


230541952573211 


X 


10" 


-2 


6.31 


X 


10" 


7 


5 


4 


4.5 





3735768429955 


1 


933923538 


X 


10" 


-5 


1 


933924400940079 


X 


10" 


-5 


4.46 


X 


10" 


7 


5 


5 


4.5 


2 


1437000387424 


7 


259841388 


X 


10" 


3 


7 


259849874157195 


X 


10" 


3 


1.17 


X 


10" 


6 


6 


5 


5.5 





8032470628900 


1 


536518877 


X 


10" 


-5 


1 


536520289551414 


X 


10" 


-5 


9.19 


X 


10" 


7 


6 


6 


5.5 


2 


5395166388946 


4 


404590776 


X 


10" 


3 


4 


404599359654937 


X 


10" 


3 


1.95 


X 


10" 


6 


7 


6 


6.5 


1 


2521347909128 


1 


148793520 


X 


10" 


5 


1 


148795883734041 


X 


10" 


-5 


2.06 


X 


10" 


6 


7 


7 


6.5 


2 


9481395894691 


2 


726943363 


X 


10" 


3 


2 


726949666682515 


X 


10" 


3 


2.31 


X 


10" 


6 


8 


7 


7.5 


1 


6997166387020 


8 


267139284 


X 


10" 


6 


8 


267168128593364 


X 


10" 


-6 


3.49 


X 


10" 


6 


8 


8 


7.5 


3 


3654788843469 


1 


713109264 


X 


10" 


3 


1 


713110726903289 


X 


10" 


3 


8.54 


X 


10" 


7 


9 


7 


8.5 





3781387034254 


1 


650193407 


X 


10" 


7 


1 


650193971262395 


X 


10" 


-7 


3.42 


X 


10" 


7 


9 


8 


8.5 


2 


1467598463840 


5 


800131555 


X 


10" 


6 


5 


800152000099428 


X 


10" 


-6 


3.53 


X 


10" 


6 


9 


9 


8.5 


3 


7890216134824 


1 


087833559 


X 


10" 


3 


1 


087830324625970 


X 


10" 


-3 


2.97 


X 


10" 


6 


10 


8 


9.5 





8445521091574 


1 


423932484 


X 


10" 


7 


1 


423939557356283 


X 


10" 


-7 


4.97 


X 


10" 


6 


10 


9 


9.5 


2 


5936449978532 


3 


998328655 


X 


10" 


6 


3 


998328388790399 


X 


10" 


-6 


6.66 


X 


10" 


8 


10 


10 


9.5 


4 


2171390613851 


6 


963970144 


X 


10" 


4 


6 


963929698290593 


X 


10" 


-4 


5.81 


X 


10" 


6 



